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An overrelaxed variant of simulated annealing is applied to the problem of maximally abelian gauge fixing. 
The superiority of this algorithm over the commonly used relaxation procedure is demonstrated. Biases on non 
gauge invariant quantities due to gauge fixing ambiguities are discussed. 



1. INTRODUCTION 

Though gauge fixing can generally be avoided in 
LGT, many physically interesting non-local quan- 
tities explicitly depend on the gauge. Promi- 
nent examples are gluon and quark propagators. 
Here, we address investigations of the dual su- 
perconductor scenario of confinement (for details 
see e.g. Ref. [1]). Maximally abelian (MA) gauge 
is believed to be most suitable for such studies. 
Fixing the MA gauge amounts to maximizing the 
functional (V = NsHes) 

with respect to local gauge transformations, 

Un,, ^ U^^ = ffn«7n,Mffi + ^, 9n G 5[/(2)/[/(l)(2) 

Besides absolute maxima - which can be degen- 
erate even beyond global gauge transformations 
and local ones within the remaining ungauged 
[/(I) subgroup - the function F(U) can have nu- 
merous local maxima. This feature resembles the 
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Gribov problem in continuum gauge theories [2]. 
Whereas degenerate global maxima, at least for 
Landau gauge, turn out to be physically equiva- 
lent [3], local maxima may lead to wrong physical 
results [4]. Therefore, in numerical practice one 
would like to find local maxima as close as pos- 
sible to the absolute ones. In this way we hope 
for smaller systematic uncertainties due to gauge 
fixing ambiguities. 

So far the relaxation plus overrelaxation (RO) 
algorithm has been employed for MA gauge fix- 
ing. In this paper we report on an implementa- 
tion of an algorithm which is successfully used in 
various optimization problems, the simulated an- 
nealing (SA) algorithm [5,6]: the functional F(U) 
is regarded as a "spin action", 

S{s) = F{U^) = ^E*^ {s^U^,,s^+,UlJ (3) 

where s„ = gliCr2,gn resemble the spin variables. 
The lattice fields Un^p, play the role of (almost) 
random local couplings. Maximizing the func- 
tional F(U^) is equivalent to decreasing the aux- 
iliary temperature T of the statistical system with 
partition function 




One starts with equilibrating this spin glass at 
high temperature. Subsequently, T is gradually 
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decreased keeping the system very close to equi- 
librium. It is evident that in the limit T — 
the system approaches its ground state, i.e. the 
maximal value of S. The advantage of SA can be 
formulated in terms of solid state physics: stan- 
dard relaxation corresponds to fast cooling and 
causes defects while the adiabatic cooling proce- 
dure of SA avoids these defects. In order to im- 
prove the movement of the spin variables through 
phase space we combine simulated annealing with 
overrelaxation (OS A). 

2. SA IMPLEMENTATION 

Our procedure consists of three steps: 1) ther- 
malization at T = 2.5; 2) gradual decrease of 
T down to T = 0.01; 3) final maximization by 
means of the RO algorithm. 

In steps 1 and 2 an overrelaxation transforma- 
tion is performed at six consecutive lattice sites 
and heatbath is applied to the seventh. Within 
step 2, every time when the heatbath update is 
applied to a site, the temperature is decreased by 
a quantum 6T. A version of this algorithm where 
sites have been visited in lexicographical order- 
ing, within subcubes of 2'* sites each, in parallel 
has also been tested successfully. The combina- 
tion of local overrelaxation and local temperature 
decreasing enables us to reduce the absolute num- 
ber of sweeps required within the cooling process 
while remaining close to an equilibrium state. 

Within the temperature range 2.5 > T > 0.1, 
6T(T) has been tuned such that the spin ac- 
tion increases about linearly with the number 
of iteration sweeps. This has been realized by 
subdividing this range into 24 intervals of width 
AT = 0.1. The corresponding differences of the 
action A5(T) = 5(T) - 5(T - AT) have been 
computed on equilibrated configurations. A5(T) 
was found to be very stable against statistical 
fluctuations among different Monte Carlo config- 
urations. An almost undetectable volume depen- 
dence and a moderate dependence on /3 was no- 
ticed. The number of sweeps (out of a fixed to- 
tal number) to be performed within each inter- 
val [T — AT, T] was computed to be proportional 
to A5(T) and, subsequently, the corresponding 
value of 6T(T) has been determined. Within the 



region 0.1 > T > 0.01, 50 additional sweeps have 
been performed. Finally, the RO algorithm has 
been applied till a convergence criterion was sat- 
isfied (about 130 sweeps in average). 

3. RESULTS AND DISCUSSION 

Our main simulations have been performed on a 
12^ lattice at /3 = 2.43. We studied the SA al- 
gorithm extensively, using three different cooling 
schedules at step 2 of our procedure: 250 (OSAl), 
500 (0SA2) and 1000 (0SA3) total sweeps. In 
addition, the standard procedure (RO) has been 
applied with the same convergence criterion. We 
collected 30 statistically independent equilibrated 
MC configurations. On each of these configura- 
tions, 10 random gauge copies have been gener- 
ated. Each of the four algorithms has been ap- 
plied to these copies. 

Let A he a gauge dependent abelian quantity. 
In the following we abbreviate A as the average 
over gauge copies and (^4) as the statistical aver- 
age. In table 1, results of the four algorithms for 
various quantities are compared with each other, 
(jj^ = [F'^—F ) denotes the variation of the maxi- 
mized value of the functional among gauge copies. 
The ideal algorithm would always give F = -Fmaxj 
i.e. asd = 0. pm is the monopole density. Wu are 
abelian i x i Wilson loops and if^^ is the abelian 
string tension. A comparison of [F) and asd ex- 
hibits the superiority of all OSA schedules over 
the standard RO algorithm, the longest sched- 
ule yielding the best results. Differences between 
the 0SA2 and 0SA3 algorithms are on the level 
of statistical errors. In means of total computer 
time spent, 0SA2 is about a factor two slower 
than RO. OSAl, still being an improvement, is 
only 30 % slower than RO. (On a CM-5, OSAl 
was found to be even slightly faster than RO.) 

It is obvious from the table that physical results 
dramatically depend on the care spent on gauge 
fixing. The functional is correlated with Wil- 
son loops and anti-correlated with the monopole 
number and abelian string tension. 

By applying state-of-the-art smearing tech- 
niques on the spatial transporters of abelian Wil- 
son loops, we have been able to compute the 
abelian ground state potential. In fig. 1 these po- 
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tentials are displayed for the RO (used by all pre- 
vious authors for MA gauge fixing) and the 0SA3 
algorithms. In addition the copy with largest 
functional among all 40 copies has been chosen 
on each configuration as our best estimate of the 
"true" maximum (circles). From this potential 
we obtain if test ~ .0478(38) as an estimate of the 
abelian string tension. (All errors have been ob- 
tained by the jackknife procedure.) From fig. 1 it 
is evident that this value will be overestimated by 
about 30% by use of the standard procedure (see 
also table 1). Even our most expensive algorithm 
0SA3 yields a value that is off by about one sta- 
tistical standard deviation (i^osAS — -0536(30)). 

The question of practical importance is how to 
obtain a string tension estimate which is reason- 
ably close to if best lowest cost. In order to 
answer this question the string tension has been 
evaluated for each algorithm using the copy with 
largest value of F out of 2, 3, etc. gauge copies. 
We came to the following conclusion: to repro- 
duce i^best 'with a systematic error of less than 
5% one can either use the 0SA3 schedule with 
2 copies per configuration or the 0SA2 sched- 
ule with 3 gauge copies which is slightly faster 
in terms of computer time. 

In order to study the scaling properties, the 
investigations have been partly repeated on a lat- 
tice with about the same physical size but smaller 
lattice resolution (16^ at /3 = 2.51). The 0SA2 
and RO algorithms have been applied to 10 copies 
on 20 configurations. The qualitative proper- 
ties are the same. Again, the string tension ob- 
tained after application of the RO algorithm is 
dramatically overestimated: K^q = .0362(15) 
versus i^best — -0305(19). Since the non abelian 
string tension at this /3 value comes out to be [7] 
K = .0337(5), this difference is relevant to the 
physical conclusions drawn. 

4. CONCLUSIONS 

It has been demonstrated that the quality of MA 
gauge fixing (in terms of the value of F(U^) 
and the scatter of results among different gauge 
copies) can be improved by applying the OSA 
algorithm even without any draw back in com- 
puter time. Our systematic study of non gauge 



invariant quantities revealed the importance of 
such an improvement and should be understood 
as a warning, not to extract the abelian potential 
and other observables without carefully investi- 
gating systematic errors, induced by the gauge 
fixing procedure. The nature of the ambigui- 
ties deserves further study. An extension of the 
method to other gauges (e.g. Landau gauge) is 
straight forward. 
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Figure 1. Abelian potentials from RO (triangles), 
OS A3 (squares) and "best" copy (circles). 



Table 1. Comparison of gauge fixing algorithms. 



RO 



OSAl 



0SA2 



0SA3 



(F) 

Pm 



0.7370(2) 

18 - 10"^ 

0.0218(2) 

0.7702(5) 

0.085(1) 

0.063(3) 



0.7383(2) 

13 - 10"^ 

0.0209(3) 

0.7716(5) 

0.090(1) 

0.057(3) 



0.7387(2) 

10 - 10"^ 

0.0207(3) 

0.7720(5) 

0.091(1) 

0.055(3) 



0.7390(2) 

8 - 10-s 

0.0204(3) 

0.7723(5) 

0.092(2) 

0.054(3) 
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